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Electron transfer organic reaction rates are considered employing the classic physical pictnre of 
Marcus wherein the heats of reaction are deposited as the energy of low frequency mechanical 
oscillations of reconfigured molecular positions. If such electron transfer chemical reaction events 
occur in the neighborhood of metallic plates, then electrodynamic interface fields must also be 
considered in addition to mechanical oscillations. Such electrodynamic interfacial electric fields in 
principle strongly effect the chemical reaction rates. The thermodynamic states of the metal are 
unchanged by the reaction which implies that metallic plates are purely catalytic chemical agents. 


PACS numbers: 31.10.+Z, 82.39.-k, 82.39.Jn, 82.75.Qt 

I. INTRODUCTION 

We consider so called electron transfer chemical reac¬ 
tions of the following sort. One starts from a molecule 
that has two parts D and A known as the donor part 
and the accepter part. A single electron passes from the 
donor D to the accepter A according to 

DA —^ D+A-. (1) 

After the electron transfer, the atomic positions within 
the molecule reorganize so as to minimize the final free 
energy. If one writes H_d as the Hamiltonian for the 
molecular configuration DA before the electron transfer 
and TIa as the Hamiltonian for the molecular configu¬ 
ration D~^A~ after the electron transfer, then the total 
model Hamiltonian may be written in the form 

"H = - (TId + T~La) + 2 — 'Ha) o’a + , 
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FIG. 1: Shown in a schematic fashion are typical diabatic en¬ 
ergy curves for Z)-states and A-states describing, respectively, 
the chemical states DA and D^A~. The adiabatic energy 
curve is the minimum of the two diabatic values. The energy 
barrier for the reaction in Eq. Q is at that reaction coordinate 
wherein the two diabatic energy curves meet. 



wherein (cti, (T 2 , cts) are the Pauli matrices describing the 
operators for the two dimensional quantum state space 
for the electron. 

The physical chemical kinetic picture of how this 
Hamiltonian works was pioneered by Marcus [THl] and 
further developed by many other workers [3-E]. The phys¬ 
ical situation[TD] is shown schematically in FIG. The 
stable free energy curves as a function of reaction coor¬ 
dinates may be defined as 

Q = Qd if Qd < Ga, 

G = Ga if Ga<Gd- (3) 

The maximum in the stable free energy takes place at 
that free energy wherein the two different free energy 
curves meet, 

G = Gma,x if Gd = G A, (4) 

at the free energy barrier upward pointing cusp. The 
chemical kinetic rate of Eq. 0 is thereby determined by 
the thermal activation flow passing over the free energy 
barrier. The minima in the free energy curves 

GD,A=^^ngD,A (5) 

yield the free energy of the reaction 

— Qjj — Qa- (6) 

On the other hand, the free energy of the barrier over 
which the electron transfer must pass is given by 

= 0max - Gd, (7) 

i.e. the kinetic reaction rate has the general thermal 
activation form 

r = 


( 8 ) 
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If the diabatic energy curves are quadratic in the reac¬ 
tion coordinates with the same force constants, then the 
resulting oscillations merely have their centers displaced 
giving rise to a reorganization potential energy A. The 
central result of Marcus for this case is the barrier energy 


^Marcus 


r(A-A0)2i 
4A 


wherein A > AQ >0. (9) 


The assumption of pure translations of a classical 
quadratic potential energy will be dropped in much of 
what follows. 

In Secjllj the chemical reaction rate will be com¬ 
puted from Eq.® on the basis of the Fermi golden rule 
r = {2Tr/h)\V\^f. The matrix element in frequency 


units is |E|/?i = 17. The density of final states gf de¬ 
termines the effective modulation time scale r = 2'Kfigf. 
Thus, r = where the condition for the validity of 
the Fermi golden rule is fir <C 1. The rigorous expres¬ 
sion for r in terms of the Hamiltonian Eq.([^ will be ex¬ 
hibited. In Sec |IIH the displaced mechanical oscillators 
will be explored and the Marcus calculation of the model 
will be discussed. In Sec|IV[ it is proposed that catalytic 
agents lower the binding energy of the product state in 
DA — D'^A~. The reason that final state binding low¬ 
ers the barrier factor and thereby increases the reaction 
rate is reviewed. This is the central mechanism we are 
proposing for how metallic catalytic agents increase the 
rate of electron transfer reactions. 

In Sec|^ we consider the reaction H 2 O —> H^OH~ 
in pure water from the viewpoint of Marcus theory. The 
metallic catalytic rate of this pure water reaction de¬ 
pends on the ordered polarized layer of water, as dis¬ 
cussed Sec |V A[ that sits on and above the metallic sur¬ 
face. In Sec |VB| the physical meaning of the experimen¬ 
tal electric dipole moment of a single water molecule is 
explored. In Sec |V C[ it is shown why the electron energy 
shift is sufficiently large to totally eliminate the energy 
barrier. 

In Sec VI the electron diabatic energy shift near the 
interface between a metal and an insulator is explored. 
The thermodynamic relations for the metal-insulator in¬ 
terface are discussed in Sec I VI Al The renormalization of 
the free energy then eliminates the barrier as explained 
in Sec |VI Bj Finally, the general reasons for the metallic 
boundaries for eliminating the barrier to electron transfer 


reactions are discussed in the concluding Sec VII 


II. REACTION RATE 


wherein there is an average over initial states and sum 
over final states. The reaction rate may thereby be writ¬ 
ten in terms of the thermal electron transfer time t. 


F = OV. (12) 

In detail 

T = 2'KTi^^piD \{fA\iD)\^ 5{EfA -Sin), 




T = 


/ oo 

Y,P^D (13) 

-00 J £ 


yielding the final electron transit time for the reaction 
DA —It is 


/ D ' 


dt, 


(14) 


wherein the average employs the thermal canonical den¬ 
sity operator pjj = for the initial states. 

If one considers a microscopic model as in Eq. ([^ , then 
Eqs.(12) and (14) are reliable for predicting chemical ki¬ 
netic electron transfer rates provided that the motional 
narrowing condition Dt <C 1 holds true. Define an inter¬ 
action V according to 


nA = nD + V, 

V{t) = 


(15) 


From Eqs.(14) and (15) one finds that 


exp ( “^ / V(s)ds 


dt, (16) 


wherein the subscript “-I-” indicates time ordering. 


III. MARCUS THEORY 


The displaced oscillator model may be described by 


— {Pk + ^kiQk — Cfc)^); 


(17) 


so that Eq.(15) reads 




9 X! ~ X! ^kikQk ■ (18) 


To compute the transition reaction rate 


The displacement reorganization energy is thereby 


r — ^DA^D+A- ! (10) 

between energy eigenstates \iD) —>■ \f A) one may employ 
the Fermi golden rule 

\{fA\iDY5{£fA-£^D), ( 11 ) 

ij 



k 


and the interaction potential Eq. reads 


V = (A - AQ) - Yj^k^kQk ■ 

k 


(19) 


( 20 ) 
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then the Marcus theory barrier free energy B for the re¬ 
action rate 


r = wherein B = 


(A-A0-W)2 
4A 


(26) 


decreases. In fact, if the binding energy W exceeds the 
upper limit 0 < -f W < A, then the barrier free en¬ 
ergy B in the reaction rate goes to zero. The situation is 
shown schematically in FIG. The renormalization in 
Eq. (24) is equivalent to switching the diabatic curve from 


Aio C lowering the reaction barrier factor B due to cat¬ 
alytic agents. Lowering the barrier free energy increases 
the reaction rate. 


FIG. 2: Shown in a schematic fashion are typical dia¬ 
batic energy curves for D-states, T-states and C-states de¬ 
scribing, respectively, the chemical states DA, A~ and 
{D'*'A~)catalytic- If the Catalytic agent increases the chemical 
binding free energy AQ, as in the dotted diabatic curve, then 
the barrier for the reaction is lowered. 


The oscillator coordinates may be considered to be classi¬ 
cal if hoJk ^ ksT in which case the equipartition theorem 
asserts that 


QkQk' — 


keT 

^k^k' 


Skk' 


( 21 ) 


wherein the average employs classical statistical mechan¬ 
ics. Thus 


V = A-A^l, 


(V-V) =J2u;tekQl = 2kBT\. 


( 22 ) 


A quasi-classical evaluation of Eq.(16) is thereby the 
Marcus result 


T = 




T = h. 


1 ^ 

<1 

1 

V XkeT 

AXksT 


r = 


(23) 


i.e. — Eqs.® and ® hold true with v = /XksT. 


V. PURE WATER REACTION 

The most simple electron transfer reaction that can 
take place in pure water is, in the notation of Eq.Q, 


H-20 = HOH 


H+OH- 


(27) 


How does a metallic surface in contact with pure water 
increase the reaction rate of Eq.( 27) required to give pure 
water the well known pH « 7 value? It appears that a 
proton is absorbed into the metal leaving an electron out¬ 
side in an ordered water layer above the metallic surface. 
That extra electron tags onto an OH bond giving rise to 
the reaction rate Eq.([27| in the form 


i70i7(liquid) - 
or equivalently 


77■''(metal) OH (liquid) 

HOH 77+ 077-. (28) 


The reaction in Eq.(28) is reversible so that the atoms 
of the catalytic metal play no direct chemical role. Fi¬ 
nally, if one inserts small spheres of metal into pure wa¬ 
ter to employ heterogeneous catalysis, then the metal 
spheres will quickly become positively charged by ab¬ 
sorbing excess protons and ejecting electrons. This is 
a well-known phenomenon, perhaps most dramatically 
displayed in the common undergraduate experiment in 
which alkali metals are dropped into water with the en¬ 
suing violent reaction which been the subject of recent 
detailed exDeriments|IIj. 


A. The Ordered Water Layer 


IV. CATALYTIC PRODUCT BINDING 

If the presence of a catalytic agent increases the bind¬ 
ing energy of the products of a chemical reaction. 


A^ ^ A^ -b W, 


wherein 


(24) 


(25) 


Since metal surfaces are in most part hydrophilic, the 
adjacent water layer is ferro-electrically ordered [12] as 
extensively studied[T2HIH| by Pollack and coworkers|16L 
I18j wherein the ordered water layer are referred to as an 
exclusion zone. What is excluded in the ordered water 
layer are positive ions and some other charged objects. 
What exists in the ordered water layer are extra electrons 
as described in the more conventional chemical symbols 
we employ above as negative ions OH~. More precisely, 
the quantitative estimates of the free energy change W 


0 < W < (A - A^l) 
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requires the physical picture of having extra electrons 
within the ordered water layer. 

The first estimate of the electron chemical potential 
renormalization of electrons within an ordered water do¬ 
main were made by Preparata employing a quantum 
electrodynamic [19] viewpoint toward ordered domain for¬ 
mation. The central physical point is as follows. If a 
single water molecule is held together by Coulomb forces 
alone, then then the energy eigenstate wave functions can 
be chosen as real and can be chosen to be eigenstates of 
the parity operator. The ground state wave function for 
a single water molecule thereby must have a zero mean 
dipole moment (0|d|0) = 0. On the other hand, the 
off diagonal matrix elements (n| d \m) with n ^ m need 
not vanish. So the physical statement that the water 
molecule has a dipole moment must refer to off diagonal 
matrix elements. If the polarization, 

= (^dj{t) 6{r-rj{t))^ , (29) 


is non-zero, then the off diagonal matrix elements of the 
dipole moment dj{t) of the molecule must appreciably 
contribute to the average in Eq.(291. Thus, in a ferroelec¬ 
tric ordered domain, the internal electronic structure is 
only partly in the ground state and partly in an excited 
state. This physical picture leads to an explanation of 
why ionic charged particles form in electrolytic water so¬ 
lutions, i.e. why it is so easy to remove an electron from 
a polarized molecule from an excited electronic state. As 
in all phase transitions, the reasons for ferroelectric or¬ 
dering reside in the collective motions of many molecules. 


B. Polarizability of a Water Molecule 

Water molecules in the ideal vapor phase have a static 
polarizability [5D] of the form 




OiT = Otoo 


SksT ’ 


aoo ~ 1.494 X 10 cm^, 
1.855 X 10"^® Gauss cm®. 


(30) 


The value of the dipole moment is often incorrectly asso¬ 
ciated with a mean dipole moment but this must thought 
through more carefully since (d) = 0. The polarizability 
is in reality 


Ut 




Pn Pm 


En — E„ 


(31) 


wherein = e^P-Eri)/kBT jg probability of being in 
state \n) and d^m = (n|d|m). From Eqs.(30) and (31) 
it follows that 




ikeT 


--E 


Pn Pm 
E — E 


(32) 


TABLE I: Work Functions into the Vacuum and into Water 


Metal 

kVvacuum 

VVwater 

Pt 

5.5 eV 

2.1 eV 

Au 

5.2 eV 

2.3 eV 

Cu 

4.7 eV 

2.1 eV 


Order of magnitude experimental determination!^ estimates of 
the work required to move an electron from the metal into the 
vacuum from the metal into and a polarized layer of water. 


wherein the restricted sum J2nm requires that energy 
differences obey \En — Em\ ^ ksT , i.e. 

= E |d„„|^ = |d|2 < (|dp) , (33) 

^—^nm 


with the inequality in virtue of the finite but small value 
of Qfoo- The root mean square fluctuation in the dipole 


moment = y |dp defines the so-called experimental 
electric dipole moment /r of a single water molecule. 


C. Electron Binding in Ferroelectric Water 

To remove an electron from an isolated water molecule 
that is initially in the ground state requires an energy 
of 1^0 ~ 12.6 eV. To excite an electron from an isolated 
water molecule that is initially in the ground state and 
finally in the first excited state requires an energy given 
by tpi « 10.2 eV. The difference (p = (po — ipi is thereby 

(j) « 2.4 eV. (34) 

As a first approximation, the work function to bring an 
electron from a metal into a ferroelectric ordered layer of 
water is Wwater ~ </'■ This work function is in qualitative 
agreement with data listed in TABLE |lj The rates of the 
pure water reaction Eq.(|27| in the neighborhood of the 
metallic surface is thereby such that the barrier would be 
virtually eliminated in virtue of A <C Wwater ~ <P- 

We note that metal catalyzed reactions in water are of 
significant current interest, especially with aim towards 
“greener” syntheses [22] . 


VI. METAL-INSULATOR CATALYSIS 

For the case of a metal catalytic agent in contact with 
an insulator, {e.g. hydrogenation of oils in the presence 
of a nickel catalyst, the Haber-Bosch synthesis of am¬ 
monia from hydrogen and nitrogen catalyzed by iron, 
or any of countless organic reactions involving nonpo¬ 
lar solvents catalyzed by metals [23|), one may apply the 
Lippmann-Gibbs eauation[24j to the interface between 
the two phases. Our sign convention is that the elec¬ 
tronic charge is negative, 

e = -|e|, (35) 
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Conductor 


hsulator 


and w denotes the charge per unit area on the charged 
sheets shown in FIG. [3] when viewed as a non-linear ca¬ 
pacitor, then the thermodynamics of such a capacitor is 
described by the Lippmann-Gibbs equation 


da = —sdT — wd^. (40) 



Dipole Moment Interface 

FIG. 3: Shown is a schematic view of the metal-insulator 
interface. The vacuum in this regard is a possible example 
of an insulator. The electrons can lower their kinetic energy 
by leaking out of the metal and into the insulator creating a 
dipole moment per unit area as pictured above. The interface 
dipole moment points from the insulator into the metal. 


The capacitance per unit area is thereby 

wherein Lt describes the effective distance between the 
charged sheets pictured in FIG. The dipole moment 
per unit area t is in a very crude order of magnitude 
given by ~ zuLt- 

B. Marcus Electron Energy Renormalization 


SO that the constant electronic chemical potential fl obeys 
/i =/r — e$ = ^-b |e|$ => ^ —/i = e$ = —|e|$ (36) 


wherein fj, is the local chemical potential and $ is the 
electrostatic potential, i.e. the electrostatic field is given 
by 

E = —grad$. (37) 

Consider the metal-insulator interface shown in FIG. |3l 
The the x-axis pointing positively normal to the interface 
from the metal into the interface, one may invoke the 
Poisson equation, 


d^$(a;) 

dx'^ 


= 47rp(a;), 


(38) 


to compute the dipole moment per unit area r employing 


vp{x)d2 


' —C30 


Attt = 


Attt = — 

d^{x) 


r . 

— oo 


dx 


dx"^ 

dx = $(oo) — $(—oo) 


47rr = -$ = ^contact < 0, (39) 


wherein Pcontact is the contact voltage across the inter¬ 
face when viewed as a possibly non-linear capacitor. The 
length scale “oo” refers to a length much larger than the 
interface dipole moment length described by Eqs.(40) 
and (411 below. 


A. Lippmann-Gibbs Equation 

If a denotes the free energy per unit area (surface ten¬ 
sion) of the interface, s denotes the entropy per unit area 


The work function to remove an electron from an in¬ 
sulating layer coating the metal may be written as 

kVinsulator “ M ~ 4” ~ 

kVinsulator — conductor | ^-hcontact |; 

kVinsulator ^ kVconductor • (42) 

The coating of a metal with an insulator lowers the work 
function. As a result, the free energy of an electron 
transfer reaction is raised via 


A^ —>■ A^ -b |eI4ontact|- 


(43) 


Since the contract potential | Pcontact | is of the order of 
volts and the reorganization energy |A/e| is of the order 
of less than than a tenth of a volt, under the renormal¬ 


ization Eq.(43) the barrier energy is eliminated. This 


allows for the explanation of the catalytic properties of 
the metal. Specific differences between different metal 
catalysts requires taking into account factors such as ab¬ 
sorption, desorption, and diffusion issues at the surfaces, 
details of the electronic structure of the metals, and sur¬ 
face preparation and is beyond the scope of this paper 
which is to simply explain why metal catalysis can be so 
spectacularly successful. We hope to return to some of 
these other matters in a future publication. 


VII. CONCLUSION 

The Marcus theory of electron transfer organic reac¬ 
tion rates has been applied to the catalytic properties 
of contact with bulk metals. The electric fields normal 
to the interface and the resulting contact potentials re¬ 
duce or even eliminate the barriers and thereby strongly 
increase the kinetic rates of reactions. The energy in¬ 
equality required for the catalytic effect, requires that 
the contact potential is large on the scale of the Marcus 
renormalization energies. 
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